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We exploit the quantum coherence between pair-produced D° and D° in i/>(3770) decays to study charm mixing 
and to make a first measurement of the relative strong phase S between D° — > K~^tt~ and D° — > i^ + 7r~. Using 
281 pb — 1 of e + e — collision data collected with the CLEO-c detector at E cm = 3.77 GeV, as well as branching 
fraction input from other experiments, we make a preliminary determination of cos<5 = 1.03 ± 0.19 ± 0.08, 
where the uncertainties are statistical and systematic, respectively. By further including other external mixing 
parameter measurements, we obtain an alternate measurement of cos 5 = 0.93±0.32±0.04, where the systematic 
uncertainty from assuming x sin <5 = has not been included. 



1. Introduction 



Recent measurements of D°-D° mixing parame- 
ters [H, H, 0, 0] highlight the need for information on 
the relative phase between the Cabibbo favored de- 
cay P° — > K~n + and the doubly Cabibbo suppressed 
decay D — » K~ir + . Here, we present a measure- 
ment that takes advantage of the correlated produc- 
tion of D° and D° mesons in e + e~ collisions. If there 
are no accompanying particles, the D°D° pair is in 
a quantum-coherent C = — 1 state. Because the ini- 



tial state (the virtual photon) has J 



PC - 1— 



there 



follows a set of selection rules for the decays of the 
D° and 5° @, i 0, B i, M EJ E M M Gl- For 
example, both D° and D° cannot decay to CP eigen- 
states with the same eigenvalue. On the other hand, 
decays to CP eigenstates of opposite eigenvalue are 
enhanced by a factor of two. More generally, final 
states that can be reached by both D° and D° (such 
as K~ir + ) are subject to similar interference effects. 
As a result, the effective D° branching fractions in this 
D°D Q system differ from those measured in isolated 
D° mesons. Moreover, using time-independent rate 
measurements, it becomes possible to probe D°-D° 
mixing as well as the relative strong phases between 
D and D° decay amplitudes to any given final state. 

In the Standard Model, D°-D° mixing is suppressed 
both by the GIM mechanism and by CKM matrix 
elements, although sizeable mixing could arise from 
new physics 1J] . Charm mixing is conventionally de- 
scribed by two small dimensionless parameters: 



y = 



J\l 2 -M l 

r 2 + ri 
r 2 + ri ' 



(i) 

(2) 



where M.\,% and I\2 are the masses and widths, re- 
spectively, of the neutral D meson CP eigenstates, 
Di (CP-odd) and D 2 (CP-even), which are defined 



as follows: 



\Di) = 
\Du) = 



\D°) + \D°) 

\D°) \D°) 
V2 ' 



(3) 
(4) 



assuming CP conservation. The mixing probability is 
then denoted by Pm = (x 2 + y 2 )/2, and the width of 
the P>° and D° flavor eigenstates is V = (T 1 + T 2 )/2. 

Many previous searches for charm mixing have fo- 
cused on D° decay times. Direct measurements of 
y come from comparing lifetimes in D° —> K + K~ 
and 7r + 7r~ decay to that in D° — > K~tt + . An in- 
direct measure of y is provided by the "wrong-sign" 
process D° — > K + tt~ , where interference between 
the doubly-Cabibbo-suppressed (DCS) amplitude and 
the mixing amplitude manifests itself in the appar- 
ent D° lifetime. These analyses are sensitive to y' = 
y cos S— x sin 5, where —5 is the phase of the amplitude 
ratio (K + 7r-\D°)/(K + 7r-\D°). Below, we also denote 
the magnitude of this ratio by r, which is measured to 
be approximately 0.06. Because S has not previously 
been measured, the separate determinations of y and 
y' above have not been directly comparable. 

In this note, we present an implementation of the 
method described in Ref. 16] for measuring y and 
cos (5 using quantum correlations at the ^(3770) res- 
onance. Our experimental technique is an extension 
of the double tagging method previously used to de- 
termine absolute hadronic Z?-meson branching frac- 
tions at CLEO-c 17]. This method combines yields 
of fully-reconstructed single tags (ST), which are in- 
dividually reconstructed D° or D° candidates, with 
yields of double tags (DT), which are events where 
both D° and D° are reconstructed, to give absolute 
branching fractions without needing to know the lu- 
minosity or D°D° production cross section. Given a 
set of input yields, efficiencies, and background esti- 
mates, a least-squares fitter [181 ] extracts the number 
of D°D° pairs produced (AT) and the branching frac- 
tions (B) of the reconstructed D° final states, while 
accounting for all statistical and systematic uncertain- 
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ties and their correlations. We employ a modified ver- 
sion of this fitter that also determines y, x 2 , r 2 , and 
r cos 5 using the following categories of reconstructed 
final states: D -> K^n ± , CP-even (5+) and CP-odd 
(S-) eigenstates, and semileptonic decays (e ± ). For 
optimal precision on 5, we also incorporate measure- 
ments of branching fractions and mixing parameters 
from other CLEO-c analyses or from external sources. 
CP violation in D and K decays are negligible second 
order effects that we ignore. 

2. Formalism 

To first order in x and y, the C-odd width 
r D ojjo(i,j) for D°D° decay to final state i/j follows 
from the anti-symmetric amplitude AAif. 

rVgo^j) oc M% = \ AiAj - AiA^ 2 

= \( l \D 2 )(j\D 1 )-( l \D 1 )(j\D 2 )\ 2 ,(5) 

where A l = (i\D°), A\ = (i\D°). The total width, 
T d oqo, is the same as for uncorrelated decay, as are 
ST rates. However, unlike the case of uncorrelated 
D°D°, we can consider the C-odd D°D° system as a 
D\D 2 pair. If only flavored final states are considered, 
as in Ref . [ItJ , then the effects of quantum correlations 
are negligible. In this analysis, we also include CP 
eigenstates, which brings additional sensitivity to y 
and 5, as demonstrated below. 

Quantum-correlated semileptonic rates probe y be- 
cause the decay width does not depend on the CP 
eigenvalue of the parent D meson, as this weak de- 
cay is only sensitive to flavor content. However, the 
total width of the parent meson does depend on its 
CP eigenvalue: = T(l =p y), so the semileptonic 
branching fraction for D\ or D 2 is modified by 1 ± y. 
If we reconstruct a semileptonic decay in the same 
event as a D 2 — > 5+ decay, then the semileptonic 
D must be a D%. Therefore, the effective quantum- 
correlated D°D° branching fractions (jF cor ) for CP- 
tagged semileptonic final states depend on y: 

T^ /e *2B s± B e (l±y). (6) 

Combined with estimates of Bi and Bs ± from ST 
yields, external sources, and flavor-tagged semilep- 
tonic yields, this equation allows y to be determined. 

Similarly, if we reconstruct a D — > K~tt + decay 
in the same event as a D 2 — ► S+, then we know the 
K~ir + was produced from a D\. The effective branch- 
ing fraction for this DT process is therefore 

FfllK. = |(5 + |P 2 )(A-^+|P 1 )| 2 
= A 2 s+ \A K -„+ + A K -^\ 2 
= A 2 s+ A 2 K ^ + \l + re^\ 2 
« B 8+ B K «(1 + Rws + 2rcos<5 + y), (7) 



where P W s = r(P>° -► K--k + )/T{D° -► R-TT+) = 
r 2 + ry' + Rm, and we have used Bs ± oc A 2 S± (1 =p y) 
and Bktt oc A 2 K _^ + (l + rycosS + rxsinS). In an 
analogous fashion, we find J r s° T /K-ir+ ~ Bs-&KirO- + 
Pws - 2r cos <5— y). When combined with knowledge of 
Bs + , y, and r, the asymmetry between these two DT 
yields gives cos<5. In the absence of quantum corre- 
lations, the effective branching fractions above would 
be B S± B K ^(1 + R WS ). 

More concretely, we evaluate Eq. with the above 
definitions of r and 6 to produce the expressions in 
Table |U In doing so, we use the fact that inclusive ST 
rates are given by the incoherent branching fractions 
since each event contains one D° and one D°. Com- 
parison of J- cor with the uncorrelated effective branch- 
ing fractions, .P unc , also given in Table Q] allows us to 
extract r 2 , rcosS, y, and x 2 . Information on Bi is 
obtained from ST yields at the i/>(3770) and from ex- 
ternal measurements using incoherently-produced D° 
mesons. These two estimates of Bi are averaged by 
the fitter to obtain .P unc . 

Table I Correlated and uncorrelated effective D°D° 
branching fractions, T C01 and T urlc , to leading order in 
x, y and r 2 , divided by Bi for ST modes i (first section) 
and BiBj for DT modes i/j (second section). Charge con- 



jugate modes are implied. 

Mode C-odd Uncorr. 

K~ir + 1 + Rws l + Rws 

S+ 2 2 

S- 2 2 

K-n + /K~n + R M Rws 

K-tt + /K + tt' (1 + Rws) 2 - 4r cos<5(rcos<5 + y) 1 + i? ws 

K~tv + /S+ 1 + Rws + 2r cos 5 + y 1 + Rws 

K~tt + /S- 1 + Rws - 2r cos 5 - y 1 + -Rws 

K~ir + /e~ 1 — ry cos <5 — rx sin 8 1 

S+/S+ 1 

S-/S- 1 

S+/S- 4 2 

S+/e- 1 + y 1 

S-/e~ 1-y 1 



3. Fit Inputs 

We analyze 281 pb -1 of e + e~ collision data pro- 
duced by the Cornell Electron Storage Ring (CESR) 
at E cm = 3.77 GeV and collected with the CLEO-c 
detector, which is described in detail elsewhere [l9[. 
We reconstruct the D° and D° final states listed 
in Table UH with 7r /?; — > 77, u> — ► 7r + 7r~7r°, and 
Kg — ► tt + h~ . Signal and background efficiencies, as 
well as crossfeed probabilities among signal modes, are 
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determined from simulated events that are processed 
in a fashion similar to data. 

Table II D final states reconstructed in this analysis. 



Type Final States 

Flavored K~ir + , K + ty~ 

S+ K+K~ , Tv+ir- , K° s -k°-k , K° L n' 
S- K° S 7v°, K° sV , K%uj 

e Inclusive Xe + u, Xe~v 



Hadronic final states without K\ mesons are fully 
reconstructed via two kinematic variables: the beam- 
constrained candidate mass, Mc 2 = \/Eq — pjjC 2 , 
where p_o is the D° candidate momentum and Eo is 
the beam energy, and AE = Eo — Eo, where Ed is 
the sum of the D° candidate daughter energies. We 
extract ST and DT yields from M distributions using 
unbinned maximum likelihood fits (ST) or by count- 
ing candidates in signal and sideband regions (DT). 

Because most K L mesons and neutrinos produced 
at CLEO-c are not detected, we only reconstruct 
modes with these particles in DTs, by demanding 
that the other D in the event be fully reconstructed. 
Ref. [2(| describes the missing mass technique used 
to identify K^tt candidates. For semileptonic de- 
cays, we use inclusive, partial reconstruction to max- 



Table III ST and DT yields, efficiencies, and their sta- 
tistical uncertainties. For DT yields, we sum groups of 
modes and provide an average efficiency for each group; 
the number of modes in each group is given in parenthe- 
ses. Modes with asterisks are not included in the standard 



and extended fits. 


Mode 




Yield 


Efficiency (%) 






25400 ± 200 


64.70 ± 0.04 


K+TT' 




25800 ± 200 


65.62 ± 0.04 


K+K~ 




4740 ± 70 


57.25 ± 0.09 


7T + 7T~ 




2100 ± 60 


72.92 ± 0.13 






2440 ± 70 


12.50 ±0.06 






7520 ± 90 


29.73 ± 0.05 






1050 ± 40 


10.34 ± 0.06 






3240 ± 60 


12.48 ± 0.04 




(2) 


4± 2 


40.2 ±2.4 


K-TT+/K+TV- 


(1) 


600 ± 25 


41.1 ± 0.2 


K^n ± /S+ (8) 




605 ± 25 


26.1 ± 0.1 


K^n ± /S- (6) 




243 ± 16 


12.3 ± 0.1 


K^/e* (2) 




2346 ± 65 


45.6 ± 0.1 


S+/S+ (9*) 




10 ±6 


12.5 ± 0.6 


S-/S- (6*) 




2 ± 2 


3.9 ±0.2 


S+/S- (12) 




242 ± 16 


7.7 ±0.1 


S+/e* (6) 




406 ± 44 


22.2 ± 0.1 


S-fe* (6) 




538 ± 40 


13.8 ±0.1 



imizc efficiency, demanding that only the electron be 
identified with a multivariate discriminant (2lT | that 
combines measurements from the tracking chambers, 
the electromagnetic calorimeter, and the ring-imaging 
Cerenkov counter. 

Table Hill gives yields and efficiencies for 8 ST modes 
and 58 DT modes, where the DT modes have been 
grouped into categories. Fifteen of the DT modes are 
forbidden by CP conservation and are not included 
in the standard fit. In general, crossfeed among sig- 
nal modes and backgrounds from other D decays are 
smaller than 1%. Modes with Kgir°ir a have approxi- 
mately 3% background, and yields for K^ir / K^tt 
and S±/S± are consistent with being entirely from 
background. 

External inputs to the fit include measurements of 
7?m, -RwSj Bk-tt+i an d Bs±i as well as an indepen- 
dent Bx^n from CLEO-c, as shown in Table Hvl The 
external i?ws is required to constrain r 2 , and thus, 
to determine cos S from rcosS. We also use the ex- 
ternal mixing parameter measurements shown in Ta- 
ble [V] The fit incorporates the full covariance matrix 
for these inputs, accounting for statistical overlap with 
the yields in this analysis. Covariance matrices for 
the fits in Ref. [13] have been provided by the CLEO, 
Belle, and BABAR collaborations. 

Table IV Averages of external measurements used in the 
standard fit. Charge-averaged D° branching fractions are 
denoted by final state. 



Parameter 


Average 


Rws 


0.00409 ± 0.00022 [22] 


Rm 


0.00017 ± 0.00039 [23J 




0.0381 ± 0.0009 [24J 


K~K + 


0.1010 ±0.0016 [25] 




0.0359 ± 0.0005 [25J 


K L TT 


0.0097 ± 0.0003 [20] 




0.0115 ± 0.0012 [24] 


K° S V 


0.00380 ± 0.00060 [24] 




0.0130 ± 0.0030 [24J 



Table V Averages of external measurements used in the 
standard and extended fits. 

Parameter Average 
y 0.00662 ±0.00211 [2,25, 26J 
x 0.00811 ±0.00334 [26] 

r 2 0.00339 ± 0.00012 [27] 

y' 0.0034 ± 0.0030 [27] 

x 12 0.00006 ± 0.00018 [2JJ 



Systematic uncertainties associated with efficiencies 
for reconstructing tracks, Kg decays, 7r° decays, and 
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for hadron identification are assigned as described in 
Ref. [29| , Other sources of efficiency uncertainty in- 
clude: AE requirements (0.5-5.5%), 7/ reconstruc- 
tion (4.0%), electron identification (1.0%), model- 
ing of particle multiplicity and detector noise (0.1— 
1.3%), simulation of initial and final state radiation 
(0.5-1.2%), and modeling of resonant substructure in 
KgTT°TT (0.7%). We also include additive uncertain- 
ties of 0.0-0.9% to account for variations of yields with 
fit function. 

These systematic uncertainties are included directly 
in the covariance matrix given to the fitter, which 
propagates them to the fit parameters. The other 
fit inputs determined in this analysis are ST and DT 
yields and efficiencies, crossfeed probabilities, back- 
ground branching fractions and efficiencies, and sta- 
tistical uncertainties on all of these measurements. 
Quantum correlations between signal and background 
modes are accounted for using assumed values of am- 
plitude ratios and strong phases that are systemati- 
cally varied and found to have negligible effect. We 
validated our analysis technique in a simulated C-odd 
D°D° sample 15 times the size of our data sample. 

4. Preliminary Fit Results 

Our standard fit excludes the 15 same-CP DT 
modes and includes the measurements in Table IIVI 
but not Table fVl In this fit, there is not enough in- 
formation to reliably determine xsin(5, so we fix it 
to zero, and the associated systematic uncertianty is 
negligible. We obtain a first measurement of cos S = 
1.03±0.19±0.08, consistent with being at the bound- 
ary of the physical region. The fit results for y, r 2 , x 2 , 
and branching fractions are consistent with previous 
measurements . 

The likelihood curve for cos<5, computed as L = 
e ~(x -x m in)/ 2 an(1 shown in Figure [TJ is slightly non- 
Gaussian. For values of | cos <5| < 1, we also show 
C as a function of | sin<5|. We integrate these curves 
within the physical region to obtain 95% confidence 
level limits of cos 5 > 0.54 and | sin 5| < 0.72. 

We also perform an extended fit that includes the 
previous measurements of y and y' in Table [Vl in ad- 
dition to all the inputs to the standard fit above. In 
this fit, we find cos 8 = 0.93 ± 0.32 ± 0.04. The sys- 
tematic uncertainty does not include the contribution 
from assuming x sin 6 = 0, which is still under study. 
From the corresponding likelihood functions shown in 
Figure [1] we determine 95% confidence level limits of 
cos S > 0.38 and |sin<5| < 0.84. 

The cos (5 uncertainty in the extended fit is larger 
than in the standard fit because of a non-linear effect. 
Most of the information on r 2 (and therefore on r) 
is provided by i?ws- Because i?ws a ls° depends on 
y ■ r cos <5, the sign of the correlation between r 2 and 
rcosS is given by the sign of y. In the standard fit, 




0.8 1.6 0.00 0.50 1.00 

cos5 |sin5| 




1 2 0.00 0.50 1.00 



cos8 |sin8| 

Figure 1: Standard (top) and extended (bottom) fit like- 
lihood functions for cos 8 and |sin<5|, including both sta- 
tistical and systematic uncertainties, excluding the uncer- 
tainty from assuming x sin S = 0. The dashed curve shows 
the Gaussian likelihood corresponding to the standard fit 
result. The hatched regions contain 95% of the areas in 
the physical regions. 



y attains a more negative central value than in the 
extended fit, where y is constrained to the precise ex- 
ternal measurements. Hence, the uncertainty on cos 5 
becomes inflated in the extended fit. 

By observing the change in {/cr 2 osS as each fit in- 
put is removed, we identify the major contributors of 
information about cos<5 to be the Kir/S± DT yields 
and the ST yields. We also find that no single in- 
put or group of inputs exerts a pull larger than three 
standard deviations on cos S or y. Moreover, removing 
all external inputs gives branching fractions consistent 
with those in Tabic HVl 

We also allow for a C-even D°D° admixture in the 
initial state, which is expected to be O(10~ 8 ) [30f , by 
including the 15 S±/S± DT yields in the fit. These 
modes limit the C-even component, which can modify 
the other yields as described in Ref. [l6| ■ In both the 
standard and extended fits, we find a C-even fraction 
consistent with zero with an uncertainty of 2.4%, and 
neither the fitted cos S values nor their uncertainties 
are shifted noticeably from the results quoted above. 



5. Summary 

Using 281 pb _1 of e + e~ collisions produced at the 
^(3770), we make a preliminary first determination of 
the strong phase 6, with cos (5 = 1.03±0.19±0.08. By 
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further including external mixing parameter measure- 
ments in our analysis, we obtain an alternate measure- 
ment of cos 8 = 0.93±0.32±0.04, where the systematic 
uncertainty from assuming x sin (5 = has not been in- 
cluded. Knowledge of 5 allows independent measure- 
ments of y and y' to be combined, thereby improving 
our overall knowledge of charm mixing parameters. 
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